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Recently it has been suggested [0 that QCD sum rule analysis can be ex- 
tended from their usual setting, which includes two and three-point functions 
[0, ^, to processes of Compton type at moderate values of the Mandelstam 
invariants s, t and u. Sum rules for a specific combination of the two invariant 
amplitudes in pion Compton scattering corresponding to different helicities, 
Hi + H2, have been given in ref. while the evaluation of relevant power 
corrections has been discussed in p. Comparison of sum rule predictions - 
in the local duality limit - with predictions from a modified perturbative cal- 
culation which includes Sudakov suppression for soft gluon exchange has 
been made in [|^ . This new approach to Compton scattering, complementary 
to the usual perturbative one, is interesting in order to investigate the tran- 
sition from non-perturbative to perturbative QCD in these processes. It also 
suggests that the non-perturbative information on this type of reactions can 
be parameterized by the lowest dimensional vacuum condensates through the 
operator product expansion of four interpolating currents. 

For the simplest Compton reaction 7r+7 —>■ 7r+7 two steps are still missing 
in completing the above programme. They are: 1) a stability analysis of the 
sum rule, and 2) the derivation of an individual sum rule for each of the 
two invariant amplitudes. Hi and H2. While here our attention is focused 
on the first point, we reserve the discussion of the second point to a future 
work. A stability analysis of the sum rule proposed in ^ is crucial in order 
to give physical justification to previous works and to allow further extension 
of the new formalism to many other similar reactions. As stated in 0, the 
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angular dependence of Compton scattering, absent in form factors, allows 
us to compare perturbative QCD and sum rule predictions in a nontrivial 
way. Our aim in this letter is first to characterize the stability region of 
sum rules for pion Compton scattering, and, second, to compare the sum 
rule predictions with the perturbative QCD approach based on the modified 
factorization formula. 

Before proceeding to study the stability of the complete sum rule, we 
recall the conclusions of IH, |^, and quote the relevant results below. A sum 
rule for the sum of the two helicities of pion Compton scattering has been 
derived by studying the following correlation function of local currents 

^af,,.x{p'l,pl, s,t) = i J d^x d'^y d'^z exp(-zpi ■ x + ip2 ■ y - iqi ■ z) 

x{0\T{vMJ,iz)MO)viix))\0), (1) 

where 

2 1 - 

Jfi = -ui^/a - -d'y^d, rj^ = uj^jad (2) 

are the electromagnetic and axial currents respectively of up and down quarks. 
These currents interpolate with the two invariant amplitudes of the scatter- 
ing process. The two photons carry on-shell momenta qi and q2, and are 
physically polarized. The momenta of the two pions are denoted as pi and 
P2, with si = pI and S2 = pi there virtualities. We also define s = {pi + gi)^, 
t = {P2~ PiY "u = {p2 — qiY for the Mandelstam invariants, which obey 
the relation s + t + u = Si + S2- The invariant amplitudes are obtained select- 
ing a specific time ordering in eq. (|1|) and projecting the two axial currents 
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onto single-pion states. Hi and H2 are isolated from the matrix element 

M,, = ij dSe-*^-^(p2|T(J^(^)J,(0)) bi) (3) 

by the expansion 

M,, = iJi(.,t)eWe(i) + H,{s,t)efe^:}\ (4) 

where e*^^-* and e*-^-* are helicity vectors defined in [|I|, 

A sum rule relates the timelike region of si and S2, where the resonant 
contribution to M^^ is located, to the so called "deep Euclidean region", 
where an operator product expansion (OPE) for the four-current correlator 
is made, through a dispersion relation. In the timelike region the spectral 
density appearing in the dispersion relation is commonly modeled by 

+Ap-*[1-0(so-pM^o-P^)] , (5) 

where the second term, the continuum contribution, which is nonvanishing 
only for p\, p\ > sq, is chosen as the perturbative spectral density A^'^''*, 
the leading term of OPE in the deep Euclidean region. Notice that the 
OPE of the spectral density includes, besides the perturbative part, non- 
perturbative power corrections proportional to the condensates of quarks and 
gluons, which are determined in close analogy with the canonical approach 
developed for the form factor case 

A Borel transform then acts on both regions in order to enhance the 
resonant contribution respect to the continuum. Given the fact that the 



dispersion relation for Compton scattering involves only a finite domain of 
the complex plane, a modified version of Borel transform 

is introduced, with C a contour of radius A^, which is kept finite in order 
to exclude the w-channel resonances from the phenomenological ansatz for 
the spectral density. In fact, can vary from sq to (s + t)/2, though it 
was set approximately to the value (s + t)/4 in ^. This modification 
introduces an extra unphysical parameter A^, in addition to the usual Borel 
mass M^, and therefore, the stability of the sum rule has to be found in the 
two-dimensional M^-A^ plane. The factor 1 — exp(— (A^ — p^)/M^) in (|^), 
which is new compared to the standard Borel transform P, ^, is to ensure 
that the transform is still regular when C crosses the branch cut in the p'f 
plane. 

The resulting sum rule for H = Hi + H2 can be expanded asymptotically 
for the large invariant 



as 



/ dsJ ds2p'"''' + —{G^) dsj dS2pS'"°M 

\J0 Jo TT Jo JO J 



where the perturbative, gluonic and quark contributions are given by, respec- 
tively, 

2560Q^V'-'(,s.Q2.,s^.,s,) 



^pert 



27s(2g2 _ Si)2(2g2 _ S2)2(4g4 _ SiS2)^(2g2s - SlS2)(4g4 - 2g2s + SiS2)2 ' 

k _ 16 {m's + 4.2 + 2M2^ + 4.^ + e) 
^ - MH ' 

with 

rP-* = (s _ g2)2(2g^5^i - gVsi + 2g"ss2 - gVs2 

-2g^siS2 - 6g^SSiS2 + 3s^SiS2) , 

+8g^°ssi + sg^s^si - losg^s^si + io4gSSi - 26g2s^si 

+8g^°ss2 + SQ^s^S2 - lOSg^s^ss + 104gVs2 - 26g2s^s2 . 

(10) 

Note that the perturbative contribution comes only from the interval (0, Sq), 
since the contribution from [sq,)^) is cancelled by that from the phenomeno- 
logical side of the sum rule. The gluon and quark condensates, (6*2) and 

((■(/'■(/')2), take the values 

^(G^) = 1.2 X lO-'GeV^ 

TT 

as{{i)il:f) = 1.8 X 10-^GeV^ . (11) 

The approximation made in the calculation of the gluonic power correc- 
tion, with the gluonic coefficient obtained by integrating its double disconti- 
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nuity along the real interval (0, A^) for each of si and S2, has been discussed 
in [0. A similar spectral representation of this coefficient has been used in 
[^] in the investigation of the form factor sum rule. In this latter case the in- 
tegration interval of such spectral representation is (0, oo), since the spectral 
density for the triangle diagram is regular for all si, S2 > and t < (see 
refs. [Q and p for details). 

Based on the formulas (0)-(0), we perform the stability analysis of 
eq. (H). The dependence of if in a wide range of the Borel mass = 4-8 
GeV^ at s = 20, \t\ = 4 and Sq = 0.7 GeV^ is shown in fig. 1. Obviously, 
all the curves located in the region marked by the vertical bars, in which 
the power corrections do not exceed 50% of the perturbative contribution, 
increase rapidly with A^. The power corrections always dominate for 
below 4 GeV^. This result indicates that there is not a stable region for H 
when the radius A^ of the Borel transform is varied. A similar behavior is 
observed for other choices of s, t and sq. As already mentioned above, such 
dependence is not present in the form factor case mainly because the vari- 
ables Si in the spectral representation for the coefficient of the gluonic power 
correction runs up to infinity [Q. 

We do not expect the strong A dependence in the sum rule for two reasons: 
1) A is an unphysical parameter, and a physical quantity like the invariant 
amplitude H should be insensitive to it; 2) the A dependence, introduced by 
the Borel transform, should cancel from both sides of the sum rule. Therefore, 
the sensitivity of our results to the radius of the Borel transform just implies 
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that the phenomenological model in (^) is too simple to maintain a stability 
of the sum rule in the M^-A^ plane. 

We propose here a modified phenomenological model in order to remove 
this spurious dependence: 



+A 



OPE 



i-e{so-pl)e{so-pl)\ , (12) 

with the continuum contribution replaced by A*-*^^, which is the same as the 
full spectral expression on the OPE side of the sum rule. This modification 
makes sense, because the region with large virtualities pf, P2 > Sq can be 
regarded as perturbative, and an OPE is allowed. With this choice we are 
requiring that the hadronic spectral density for the continuum from sq to 
truncates not only the perturbative part (the lowest order contribution), but 
also the power corrections, of the OPE side of the sum rule. Since the contri- 
butions from the region (sq, A^) have been removed by the above cancellation, 
Sj's never reach the upper bound A^, and the remaining A dependent factors 
1— exp[— (A^ — Si)/M^] from the transform (^) can be dropped. Therefore, the 
overall dependence on the radius of the Borel transform disappear completely 
from the sum rule. 

Using eq. (|12|), (||) is modified to 
's{s + t)\ 



-t 



Jo Jo n Jo Jo J 

+C'^''^'^7fa,{{4j4jf) , (13) 



with the modified quark contribution 



^quark 



.(14) 



9 MH^ 
Note the change of the upper bound from to sq in the integral for the 
gluonic power correction, which is due to the cancellation from the phe- 
nomeno logical side. 

Before studying the new sum rule (|I3D, we shall examine how the modified 
parametrization for the continuum affects the sum rule calculation of the pion 
form factor. This is a significant check for the process we aim to discuss, since 
it has been shown |I| that Compton scattering has a strong similarity to the 
corresponding form factor case at moderate s and t and at a fixed angle. The 
modified sum rule for pion form factor can be derived easily based on [^] : 



JO Jo TT Jo Jo J 

, (15) 



with 



pert 



gluon 



3g" 1 

167r2 57/2 
1 



48M4 



35(Si + S2 + g')(Si +S2 + 2g2) _ 52 _ 5g2^^^ ^ ^ q2^ 

Q\si + S2) + {si-S2f-25-2Q'' 



p/2 



^ AQ\si + S2 + Q^) _^si + S2 + Q^ 
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M251/2 



^quark 

6 



52 



81M4 \ 13M2 J ' 

(Si + S2 - t)^ - AS1S2 



(16) 



Note that the expression for the quark power correction C!,^^^^^ based on the 
modified phenomenological model is the same as in p, Q . A Borel transform 
picks up only the residue of the pole l/PiP^ in the quark power correction, 
and there is not such a pole term in the integral from sq to on the phe- 
nomenological side of the sum rule. Therefore, the quark part in A*-*^^ does 
not contribute, when the Borel transform is applied. There is not cancellation 
from the phenomenological side for the quark power correction. The differ- 



ence of the modified expression ( |T4|) in the case of pion Compton scattering 
from the corresponding one in (|^) is due to the neglect of the suppressing 
factors mentioned above, which results in the change of the spectral density. 

Results for F^^ obtained from the stability analysis of (|T5|) are shown 
in fig. 2 with the best values of sq = 0.7 and = 2 GeV^ substituted 
into the sum rule. It is evident that the behavior of in is indeed in 
good agreement with experimental data [|| even after these amendments, 
and consistent with those derived in P, which are based on the standard 
phenomenological model (H). 

Now we proceed to the stability analysis of (|TB|) following a method similar 
to 0, ^. Since the A dependence has been removed completely, we concen- 
trate simply on the variation of H with respect to IlP. The dependence 
of H for So = 0.5 - 0.7 GeV^ at s = 20 and \t\ = 4 GeV^ is displayed in 
fig. 3. The region on the right-hand side of vertical bars is the one where 
the power corrections do not exceed 50% of the perturbative contribution. It 
is obvious that as Sq = 0.6 GeV^ there is the largest interval, in which 



10 



H is approximately constant. Therefore, sq = 0.6 GeV^ is the best choice 
which makes both sides of the sum rule most coincident. This value of the 
duality interval is close to that given in the form factor case, and consistent 
with its conjectured value 0.7 GeV^ in 0. Different sets of s and t have 
been investigated. The best value of sq does not vary significantly, and H is 
almost constant within the range 2 < < 6 GeV^. 

Results for H at different scattering angles of the photon, 9, sin(^^/2) = 
— t/s, with So = 0.6 and = 4 GeV^ are exhibited in fig. 4, where \H\ 
denotes the magnitude of H . Basically, they show a similar dependence 
on angles and momentum transfers \t\ to those derived using local duality 
approximation |^, but with the magnitude lower by 20% only at small \t\. 
These predictions are compared to the perturbative predictions obtained from 
the modified factorization formula [0]. The transition to perturbative QCD 
at about \t\ = 4 GeV^ and —t/s = 0.5 {6 = 40°), where the perturbative 
contributions begin to dominate, is observed. Sum rule results are always 
smaller for —t/s = 0.6 {6 = 50°), and always larger for —t/s = 0.2 {6 = 15°), 
than the perturbative results. This is also consistent with the conclusion in 

i- 

We have analyzed in more detail the sum rule which describes the be- 
haviour of the sum of the two helicities of pion Compton scattering close 
to the resonant region, and its dependence on the two parameters, and 
A^, which characterize the modified Borel transform. We have seen that the 
strong dependence of the sum rule on the radius of the transform can be re- 
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moved consistently under the assumption that resonances from mass higher 
than the duahty interval sq contribute equally both to the phenomenological 
side and to the OPE side. The results for the two helicities are found to be 
stable under the new phenomenological parametrization for the continuum 
in a wide variation of the Borel mass, and is characterized by a local duality 
interval which takes a value similar to the form factor case. A more detailed 
discussion of these issues on individual sum rules for the two helicities will 
be considered elsewhere 0. 
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Figure Captions 



Fig. 1 Dependence of if on at s = 20, |t| = 4 and sq = 0.7 GeV^ for 
= 4 (solid line), 6 (dashed line), and 8 GeV^ (dotted line). 

Fig. 2 Dependence of on derived from the modified phenomenological 
model dll) with sq = 0.7 and = 2 GeV^ (solid line). Resuhs from ref. g 
(dashed line) and from ref. (dotted line), and experimental data (dots) 
are also shown. 

Fig. 3 Dependence of if on at s = 20 and \t\ = 4 GeV^ for (a) sq = 0.7 
GeV2, (b) So = 0.6 GeV2, and (c) sq = 0.5 GeV^. 

Fig. 4 Dependence of |t| |ii| on \t\ derived from the full analysis of QCD sum 
rules with sq = 0.6 and = 4 GeV^ (sohd lines) for (a) -t/s = 0.6 (9 = 
50°), (b) -t/s = 0.5 {9 = 40°), and (c) -t/s = 0.2 (9 = 15°). Corresponding 
results from the modified perturbative QCD calculation (dashed lines) are 
also shown. 
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